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Abstract — This paper studies the achievable rate region of the 
two-user single-input-single-output (SISO) Gaussian interference 
channel, when the improper Gaussian signaling is applied. Under 
the assumption that the interference is treated as additive 
Gaussian noise, we show that the user's achievable rate can 
be expressed as a summation of the rate achievable by the 
conventional proper Gaussian signaling, which depends on the 
users' input covariances only, and an additional term, which is a 
function of both the users' covariances and psendo-covariances. 
The additional degree of freedom given by the pseudo-covariance, 
which is conventionally set to be zero for the case of proper 
Gaussian signaling, provides an opportunity to improve the 
achievable rate by employing the improper Gaussian signaling. 
Since finding the optimal solution for the joint covariance 
and pseudo-covariance optimization is difficult, we propose a 
sub-optimal but efficient algorithm by separately optimizing 
these two sets of parameters. Numerical results show that the 
proposed algorithm provides a close-to-optimal performance as 
compared to the exhaustive search method, and significantly 
outperforms the optimal proper Gaussian signaling and other 
existing improper Gaussian signaling schemes. 

I. Introduction 

The capacity of the two-user Gaussian interference channel 
(IC) has been an open problem for a long time [ 1 1 . Recently, 
a significant progress has been made in [2], where it is proved 
that a particular Han-Kobayashi type scheme can achieve 
within one bit to the information-theoretical capacity. A key 
technique in the Han-Kobayashi scheme is to split each user's 
transmit signal into a common message, which is decodable 
at both receivers, and a private message, which is decodable 
at the intended receiver only. Since this capacity-approaching 
technique requires signal-level encoding/decoding coopera- 
tions among the users, which are difficult to implement, a 
more pragmatic approach is to implement single-user detection 
at the receivers by treating the interference as noise |3)-|]5). 

There has been a great deal of research on characterizing 
the Pareto boundary of the achievable rate region for the 
Gaussian IC with the interference treated as noise BJ-JT)- 
A Pareto boundary consists of all the rate-tuples at each of 
which it is impossible to improve one particular user's rate, 
without simultaneously decreasing the rate of at least one of 
the other users. A common approach for such characterizations 
is via solving a sequence of weighted sum-rate maximization 



(WSRMax) problems (5). However, as pointed out in [4], the 
WSRMax approach cannot guarantee the finding of all Pareto- 
boundary points due to the non-convexity of the achievable 
rate set. An alternative method based on the concept of rate 
profile was proposed in BJ, which is able to characterize 
the complete Pareto boundary for the multiple-input-single- 
output IC (MISO-IC) . Furthermore, the rate-profile approach 
generally results in optimization problems that are easier to 
handle than the conventional WSRMax problems HI. 

However, all the aforementioned works are restricted to 
proper Gaussian input signals, for which the second-order 
statistic is completely specified by the covariance matrix 
(under the zero-mean assumption). On the other hand, for the 
more general improper Gaussian signaling, an extra parameter 
called pseudo-covariance is required for the complete second- 
order characterization of the complex-valued input signals 
1 8 1 — [ 1 1 . This extra parameter provides a new opportunity to 
improve the achieve rates of the Gaussian ICs. For instance, 
it was shown in fTT) that improper signaling is beneficial 
in improving the degrees-of-freedoms (DoF) performance for 
the three-user single-input-single-output IC (SISO-IC) with 
time-invariant channel coefficients. Since the DoF metric is 
meaningful only at the asymptotically high signal-to-noise 
ratio (SNR), in this paper, we are interested in characterizing 
the achievable rate region with improper Gaussian signaling 
at any finite SNR value. For the purpose of exposition, we 
consider the simple two-user SISO-IC in this paper. 

The achievable rate region for the two-user SISO-IC with 
improper Gaussian signaling has been studied in 1 12|13) , based 
on the equivalent 2x2 real-valued multiple-input-multiple- 
output (MIMO) channel obtained by separating the real and 
imaginary parts of the channel coefficients. In [12], a rank-1 
signaling scheme was proposed, i.e., the transmit covariance 
matrices in the equivalent MIMO channel are restricted to 
be rank one. This is equivalent to transmit purely real or 
purely imaginary signals over the original complex-valued 
SISO channel. In fl3| , the rate region was obtained based 
on an exhaustive search over the input covariance matrices. 
It was shown that by exploiting the symmetry property of 
the covariance matrices, the search space can be reduced to 
the 4-dimension; however, the resulted search is still of high 



complexity. In contrast to the above prior works, in this paper, 
we adopt the complex-valued SISO channel model to gain 
some new insights. Based on existing results on improper 
complex random vectors (RVs), we show that the achievable 
rate with improper Gaussian signaling for the two-user SISO- 
IC can be expressed as a summation of the rate achievable by 
the conventional proper Gaussian signaling, which depends 
on the users' input covariances only, and an additional term, 
which is a function of both the users' covariances and pseudo- 
covariances. By applying the rate-profile technique, we then 
propose an efficient algorithm to optimize the covariances and 
pseudo-covariances to enlarge the achievable rate region. 

The rest of this paper is organized as follows. Section [II] in- 
troduces the system model and some preliminaries. Section III 
presents the problem formulation and the proposed algorithm. 



Numerical results are provided in Section IV Finally, we 
conclude the paper in Section [V] 

II. System Model 

Consider a two-user SISO-IC, where each transmitter is in- 
tended to send one independent message to the corresponding 
receiver. The input-output relationship is given by 

Vi = huXi + h 12 x 2 + (1) 
V2 = h 2 \X\ + h 22 x 2 + n 2 , (2) 

where y\ and y 2 are the received signals at receiver 1 and 
2, respectively; h rt = \h rt \e^ rt 1 r 1 1 = 1,2 denotes the 
complex channel coefficient from transmitter t to receiver r 
and (j) r t represents its phase; m and n 2 are the zero mean 
circularly symmetric complex Gaussian (CSCG) noises with 
identical variance a 2 , denoted by nx,n 2 ~ CAf(0, cr 2 ); and 
x\ and x 2 are the independent signals from transmitter 1 and 
2, respectively. Different from the conventional setup where 
proper Gaussian signaling is assumed, in this paper, x\ and 
x 2 are zero-mean complex Gaussian random variables which 
can be improper. 

A. Preliminary: Improper Random Vectors 

For a zero-mean RV z 6 C™, the covariance matrix C z and 
the pseudo-covariance matrix C z are defined as [81 

C z ^E(zz H ), C z ^E(zz T ), (3) 

where (-) T and (-) H represent the transpose and complex- 
conjugate transpose, respectively. 

Definition 1. jj^V: A complex RV z is called proper if 
its pseudo-covariance matrix C z vanishes to a zero matrix; 
otherwise it is called improper. 

Define C z as the covariance matrix of the augmented vector 
[z T z* T ] T , where (■)* represents the complex conjugate 
operation, i.e., 
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(4) 



Theorem 1. [9J: C z and C z are a valid pair of covariance 
and pseudo-covariance matrices, i.e., there exists a RV z with 



covariance and pseudo-covariance matrices given by C z and 
C z , respectively, if and only if the augmented covariance 
matrix C z is positive semidefinite. 

Theorem 2. [9]: The differential entropy of a complex 
Gaussian RV z with augmented covariance matrix C z is given 
by 



Mz) = -log[(7re) 2rl detC z ]. 
B. Achievable Rate with Improper Gaussian Signaling 



(5) 



For a scalar complex random variable z, we use C z and C z 
to denote the covariance and pseudo-covariance, respectively. 
Then for the zero-mean input Gaussian signals x\ and x 2 , we 
have 



C Xl = E(xixl), C Xl =E(xia;i) 
C X2 = E(x 2 x 2 ), C X2 = E(x 2 x 2 ) 



(6) 
(7) 



Note that C Xl and C X2 are nonnegative real numbers equal to 
the power values of the transmitted signals, while C Xl and C X2 
are complex numbers in general. From Theorem [T] it is easy 
to verify that the following conditions are both necessary and 
sufficient for C Xl and C Xl (or C X2 and C X2 ) to be a valid pair 
of covariance and pseudo-covariance for a random variable X\ 
(or x 2 ) 



\C Xl \ < C Xl , \C X2 \ < C X2 



(8) 



Next, we derive the rate expression in terms of the input 
covariances and pseudo-covariances. With Gaussian input X\ 
and x 2 , it is well known that the received signals y r , r = 1,2, 
are also Gaussian. The covariance and pseudo-covariance of 
y r can be obtained as 



C Vr = E(y r y* r ) = \h rl \ 2 C Xl + \h r2 \ 2 C X2 + a 2 , 

C Vr = E (yrVr) = hl 1 C Xl + h 2 r2 C X2 . 



(9) 



With the conditions given in (fSD, it can be verified that 
|Cj, r | 2 < C 2 r is satisfied. Then from Theorem |ij the dif- 



ferential entropy of y r is given by 



KVr) = ^log(Tre) 2 det 



Cy r Cy 

cz c v 



\og(neC yr ) + -log(l-C; 2 \C yr \ 2 ). 

n r , which is the interference plus noise 
term at receiver r, where r = mod (r + 1,2). Then 



C s 



~C X - + a , C Sr 



' b rr x ~r' 



(10) 



Similarly, the differential entropy of s r can be obtained as 

1, 



h(s r ) = log(7reC s 



-\og{\ ~C~ 2 \C S 



(11) 



Under the assumption that interference is treated as additive 
Gaussian noise and perfect channel knowledge is known at 



all terminals, the achievable rate at receiver r with improper 
Gaussian signaling can be obtained as 



R r = I(x r ; y r ) = h(y r ) - h(s r ) 

i i \h rr \ 2 C x s 1 

v v ' 

flr pc '(c xl ,c X2 ) 



log 



i-c-'\c Vr \* 

1 - C7 2 \C S \ 2 



(12) 

Equation ( |T2] > clearly shows that, compared to the conventional 
proper Gaussian signaling, the achievable rate with improper 
Gaussian signaling has an additional term, which is a function 
of both the covariances and pseudo-covariances of the input 
signals. By setting C Xl and C X2 both to be 0, ( fT2] > reduces 
to the well-known rate expression for the proper Gaussian 
signaling. 

III. Achievable Rate Region with Improper 
Gaussian Signaling 

The achievable rate region for the two-user IC is defined to 
be the set of rate-pairs for both users that can be simultane- 
ously achieved under a given set of transmit power constraints 
for each transmitter, denoted by Pi, P 2 , i.e.,: 



by jointly optimizing the covariances and pseudo-covariances 
is quite involved. We therefore propose a suboptimal solution 
in this paper with separate optimization of the covariances and 
pseudo-covariances. Specifically, with pseudo-covariances set 
to zeros, (PI) is shown to reduce to a convex optimization 
problem, which can be efficiently solved with linear program- 
ming (LP). On the other hand, for fixed covariances, (PI) is 
shown to be equivalent to solving a finite number of second- 
order cone programming (SOCP) problems, from which the 
pseudo-covariances can be optimally solved efficiently. 

A. Covariance Optimization 

When restricted to proper Gaussian signaling with C Xl = 
and C X2 — 0, (PI) reduces to 

(Pl-a): max. r 



s.t. log(l 



log(l 



\hn\ 2 C Xl 
r 2 + \h 12 \ 2 C x 



'•22 



"C x 



-) > ar, 



) > (l-a)r, 



(13) 



(14) 



U 



{(r l! r 2 ):0<r 1 < R u < r 2 < R 2 }, 



C^<Pi,\C xl \<C xl 
C X2 <P 2 ,\C X2 \<C X2 

where R\ and R2 are given by ( [12) . 

Definition 2. \7y: A rate-pair (ri,r 2 ) is Pareto optimal if 
there is no other rate-pair (r' 1; r' 2 ) with (r[, r 2 ) > (7*1, r 2 ) and 
( r ii r ' 2 ) 7^ ( r i> r 2), where the inequality is component-wise. 

For the achievable rate region 1Z with improper Gaus- 
sian signaling, we adopt the rate-profile technique in |4j to 
characterize the Pareto optimal rate-pairs. Specifically, any 
Pareto optimal rate-pair of 1Z can be obtained by solving the 
following optimization problem with a particular rate profile 
denoted by (a, 1 — a): 

(PI): max. _ R 

jC^e 2 >C' X 2 '"^ 

s.t. Rr a (C Xl ,C X2 ) + ^log 1 ~^j < H > aR, 



1. 



i?r per (c, 1 ,c7 ;C2 ) + -iog 



'i-c^\c si \ 2 

1 — Cy 2 \Cy 2 I 



1 



Cs 2 2 \C S2 1 2 



> (1 - a)R, 



0<C X1 < P u 0<C X2 < P 2 , 



a 2 + \h 21 \ 2 C Xl 
0<C Xl <P l5 0<C X2 <P 2 . 

(Pl-a) is non-convex and hence cannot be solved directly. 
However, for any fixed value r, (Pl-a) can be transformed 
to the following LP feasibility problem: 

(Pl-a'): Find C Xl ,C X2 

s.t. \h n \ 2 C Xl > (a 2 + \h 12 \ 2 C X2 )(e ar - 1), 

\h 22 \ 2 C X2 > (a 2 + \h 21 \ 2 C Xl )(e^ r - 1), 
< C X1 < P u 0<C X2 < P 2 . 

(Pl-a') can be efficiently solved with existing algorithms such 
as the simplex method [14) . If r is feasible to (Pl-a'), then 
it follows that the optimal solution to (Pl-a) satisfies r* > r; 
otherwise, r* < r. Thus, (Pl-a) can be efficiently solved by 
solving (Pl-a') with different values for r, together with the 
bisection method for updating r [15]. Denoting the optimal 
solution to (Pl-a) as {r*, C* , C* 2 }, then it can be verified 
that the constraints ( fl3] i and ([14} will be both active, i.e., 

<° per (C- ,C-) = ar*,^^,^) = (1 -a)r\ 

(15) 

B. Pseudo-Covariance Optimization 

In this subsection, (PI) is optimized over the pseudo- 
covariances C Xl and C X2 , by fixing the covariances to C£ 
and C* 2 obtained by solving (Pl-a). The resulted problem is 
formulated as 



I a 



— CjCl! 1^2 I 



< c 2 



where a € [0,1] denotes the target ratio between user l's 
achievable rate and the users' sum-rate, R. Denote the optimal 
solution to (PI) as R*, then the rate-pair (aR*, (1 — a)R*) 
must be on the Pareto boundary corresponding to the rate 
profile given by a. Thereby, by solving (PI) with different a 
values between and 1, the complete Pareto boundary for the 
achievable rate region 7Z can be found. However, solving (PI) 



(Pl-b): _ max. R 

C X1 ,C X2 ,R 



s.t. RT^iC^C^ + ^o; 

RT ve \c* xl ,c* X2 )- 11,,. 



2 S 1-C S - 2 |C S1 P 



> aR, 



1 C S2 |C S2 | 2 



\c xi \ 2 <c* x 2 , \c X2 \ 2 <c* 
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Fig. 1: Achievable rate region for SNR = dB 



where C yi ,C Sl , G V2 and C S2 are the corresponding covariance 
terms with input covariances C* and C* 2 . Again, if a given 
R is achievable for certain C Xl , C X2 , then the optimal solution 
to (Pl-b) satisfies R* > R; otherwise, R* < R. This enables 
solving (Pl-b) via solving a set of feasibility problems, each 
for a fixed value R. Moreover, R can be updated with a simple 
bisection search fl5|. Substituting ( fT3] l into (Pl-b), it can be 
easily obtained that {C Xl = 0, C Xa = 0, R = r*} is feasible to 
(Pl-b). Therefore, R* > r* is satisfied, i.e., with our proposed 
separate covariance and pseudo-covariance optimization, the 
sum-rate corresponding to the rate profile given by a with 
improper Gaussian signaling is no smaller than that obtained 
with the optimal proper Gaussian signaling. 

Next, we present the algorithm for the feasibility problem 
resulting from (Pl-b) for a given target R. Substituting |[9), 
(jTOf and ( fT3j ) into (Pl-b) and after some simple manipulations, 
the feasibility problem for a given R can be formulated as 



(Pl-b') _min. 

Cxi fCm 







s.t. a^hltC^ +h 2 12 C : '- 



a 2 \h 21 C Xl 

\c xi \ 2 <c%, 



%2 I 

^22^X2 I 



bi < \C J2 \ 2 , 
b2<\C Xl \\ 



\c x 



where a\ = 



< C* 



(16) 
(17) 
(18) 
(19) 
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i«2ir 

2(l-a)(.R-r*) 



Since R* > r*, we can assume that R > r* without loss of 
optimality. Then it follows that p± > 1,(32 > 1, &i > and 
b 2 > 0. 

(Pl-b') is non-convex and hence cannot be solved directly. 
Next, we show that it can be efficiently solved via solv- 
ing a finite number of SOCP problems. First, observe that 



an arbitrary common phase rotation can be added to both 
pseudo-covariances without affecting the feasibility of (Pl- 
b'). That is, if {C Xl ,C X2 } is feasible for (Pl-b'), then so is 
{C Xl e^, C X2 e JW }. Therefore, without loss of generality, we 
may choose ui so that C X1 is real and nonnegative. Denote 
the magnitude and phase of C X2 by t and 6, respectively, i.e., 
C X2 = t&* b . Then for any fixed value of 9, (Pl-b') can be 
transformed into a SOCP feasibility problem given by 





(Pl-b"): min. 

s.t. 



MhlMC Xl } + h 2 12 te^ 9 ) 



< t, 



<^{C K1 }, 



%{C X1 } = 0, n{c Xl }<c* Xi , t<c* r 

(Pl-b") is convex and hence can be efficiently solved with the 
standard interior point algorithm fT3) , or by existing software 
tools such as CVX (16). 

Theorem 3. The feasibility problem (Pl-b') can be optimally 
solved by solving a finite number of SOCP problems (Pl-b"), 
each for a fixed value 9, where 9 can be restricted to the 
following discrete set: 



9 € {7 



),7T + 2( 



22 )} ue A u &t 



where 0_4 and Og are the solution sets for 9 to the following 
equations: 



9t 



h\ 2 te> 



a x \h\ x C% x 
a 2 \h 2 21 C* Xi + h* 22 te 



0)2 



012 



b 2 



*2 



a 1 \h\ 1 X 1 



hf 2 C^e 



je\2 



bi = Cf 2 



a2\h 2 21 X 1 +h 2 22 C*e^\ 2 + b2=X 2 1 



(20) 



(21) 



Proof: Please refer to Appendix [A] ■ 
Theorem [3] can be intuitively interpreted as follows. For 
the feasibility problem (Pl-b'), if the constraint ( |T~6] > is more 
"restrictive" than ( fT7] >, then 9 should have a value such that 
the left hand side of ( fTo} is minimized. This corresponds to 
9 = 7r+2(</)ii— <fii2) so that h 2 1 C Xl and h\ 2 C X2 are antiphase. 
Similar interpretation for 9 = tt + 2(4>2i — <\>22) can be made. If 
both ( fT6| ) and ( fTTj ) are equally "restrictive", a feasible solution 
tends to make both constraints satisfied with equality, as given 
by ( |20] > and ( |2"T| i. The elements in 6^ and 9g can be easily 
obtained as shown in Appendix [B| 

IV. Numerical Results 

In Fig. [T] and Fig. |2j the achievable rate regions (prior 
to any time-sharing of achievable rate-pairs) of the proposed 
improper Gaussian signaling scheme at SNR=0 versus 10 dB 
are compared with other schemes, including the optimal proper 



Gaussian signaling scheme presented in Section III-A the 
optimal improper Gaussian signaling obtained with the ex- 
haustive search method [131, and the rank-1 scheme with both 



Improper, exhaustive 
* Improper, rank-1 with MMSE 
o Improper, rank-1 with ZF 
-e- - Proper 
-A — Improper, proposed 




by H 



The 
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Fig. 2: Achievable rate region for SNR = 10 dB 



zero-forcing (ZF) and minimum-mean-square error (MMSE) 
beamforming [12|. The channel matrix for both plots is given 
r 1.5718 - 1.2863? -1.2984 + 0.7032il 
.-0.2847 + 0.6700i 0.7802 - 0.615H . 
(r, £)fh element of H is h rt , which is the channel coefficient 
from transmitter t to receiver r. Both figures reveal that the 
achievable rate regions have been significantly enlarged with 
improper Gaussian signaling. It is also observed that our pro- 
posed algorithm, albeit being sub-optimal due to the separate 
optimization of covariances and pseudo-covariances, performs 
quite close to the optimal improper Gaussian signaling by the 
exhaustive search. Another interesting observation is that for 
this particular channel realization, the Pareto boundary points 
of the achievable rate region with time-sharing (by taking 
the convex-hull operation over all the achievable rate-pairs 
without time-sharing) with improper Gaussian signaling can 
be obtained by the time-sharing between the two single-user 
maximum rate points, and the north-east rate corner point of 
the rank-1 scheme with MMSE beamforming. 

V. Conclusion 

This paper studies the achievable rate region of the two-user 
SISO Gaussian IC, when the improper Gaussian signaling is 
applied. It is shown that the achievable rate can be expressed 
as a summation of the conventional rate expression by proper 
Gaussian signaling, and an additional term that is a function 
of both the input covariances and pseudo-covariances. An effi- 
cient algorithm is proposed to obtain enlarged achievable rate 
regions by optimizing the covariances and pseudo-covariances 
separately with improper Gaussian signalling. 

Appendix A 
Proof of Theorem[3] 

For notational convenience, in the sequel, we use X\ and 
X 2 to denote C Xl and C X2 , respectively. First, the following 
proposition shows that to solve (Pl-b'), we may consider 
exterior solutions only, i.e., the solutions at which at least one 
of the inequality constraints is active. 



Proposition 1. If{Xi,X 2 } is feasible to (Pl-b') with \X X \ < 
C* and \X 2 \ < C* , then there exists another feasible 
solution {X[,X 2 } with \X[\ = C* x or \X' 2 \ = C* 2 . 



Proof: Let t = minlp^j, r^Jr}- Then r > 1. Define 
X[ = tXi,X 2 = tX 2 . Then it is easy to verify that the 
constraints in ( fT8) > and ( fT9] > are satisfied, i.e., \X[\ < C* x and 
\X 2 1 < C* 2 . Furthermore, at least one of them are satisfied 
with equality. The constraint in ( fT6] > is also satisfied since 

a^h^X'i + h\ 2 X' 2 \ 2 +b x = tVKi^i + h\ 2 X 2 \ 2 + b x 



(a) (b) 

< T 2 ( ai |^ll*l+^2*2| 2 +M < T 2 \X 2 \ 2 = \X' 2 \ 2 , 



where (a) is satisfied since r > 1 and b\ > 0, (b) is true 
since {X X ,X 2 } is feasible to (Pl-b'). Similarly, ( fT7| i is also 
satisfied. Therefore, {XJ,^} is a feasible solution to (Pl-b') 
with at least one of the inequality constraints being active. ■ 
Next, we derive Theorem [3] using the KKT conditions, 
which are necessary optimality conditions for the constrained 
optimization problem (Pl-b') p5| . For notational convenience, 
denote the inequality constraints (p~6]>— (|T9]> by fi < 0, f 2 < 
0,hi < and h 2 < 0, respectively. Denote Ai, A2, /ii, [i 2 
as the corresponding dual variables, respectively. Then the 
Lagrangian function of (Pl-b') is given by 

L{X 1 ,X 2 ,XxA2,Hi,^) = Xi{a 1 (\h 11 \ 4 X* 1 X 1 + \h 12 \ i X* 2 X 2 
+ 2di{h 2 * 1 h 2 2 X* 1 X 2 }) + 61 - X|X 2 } + \ 2 {a 2 (\h 21 \ i X*X 1 



+ \h 22 \ A X*X 2 + 2ft{h 2 2 lh 2 22 X*X 2 }) + b 

+ MiC^Xi - cfj + » 2 (x*x 2 - c: 2 ). 



X 



(22) 



For {X x , X 2 } to be a solution to (Pl-b'), the following KKT 
conditions must be satisfied [15|: 

1) Dual feasibility: A x > 0, A 2 > 0, m > 0, fi 2 > 0. 

2) Zero derivative: The derivatives of the Lagrangian func- 
tion (p2|i with respect to the primal variables are zero: 



dL 

dX* 



= 



X 2 [(ai|/ii2| 4 -l)Ai + A 2 a 2 |/i 2 



M2 



dL 

dX* 



= Xi(\i axh^h^ +\ 2 a 2 h 22 h 2l ) 

Vi Vi 

= =>X X [{a 2 \h 21 \ 4 - 1)A 2 + A iai |/i n | 4 + /n] 

Cl 

= -X 2 (Ai aih 2 \h\ 2 +A 2 a 2 h\\h\ 2 ), 
V* v 2 * 

-X 2 c 2 =X 1 (\ 1 V 1 + \ 2 V 2 ), 
X x ci = -X 2 (A 1 F 1 * + A 2 F 2 *). 



(23) 
(24) 



3) Complementary slackness: A1/1 = 0, A 2 / 2 = 0,//i/ii 
0,/j, 2 h 2 = 0. 



As discussed previously, without loss of generality, Xi can be 
assumed to be a nonnegative real number. The following cases 
are then considered to derive the possible phases of X 2 : 
• Case I: /1 = and f 2 7^ 0. Then from the complementary 

slackness condition, Ai > and A2 = 0. Substituting 

them into ( |24| >, we have 

Xi(Aiai|/in| 4 + Hi) 

Since A x > 0, fii > 0, ai > and Xi > 0, the phase 
9 of X2 equals to that of V\ rotated by 7r, which is 7r + 
2(0n — $12) since V\ = aihf 2 h 2 1 . 

. Case II: fi^O and f 2 = 0. Then Ai = and A 2 > 0. 
Similarly, by using ( |23j ), we have 8 = it + 2(02i — 022)- 

. Case III: f 1 ^ and f 2 ^ 0, then Ai = and 
A2 = 0. By substituting them into |23) and p4| , we have 
M2AT2 = and [i,\X\ = 0. This means either X2 = 0, 
Xi = 0, or /ii — 0, /i2 = 0. The former case is trivially 
corresponding to proper Gaussian signaling. Furthermore, 
Proposition [T] suggests that we may consider the exterior 
solutions only, i.e., either hi = or h% = is satisfied. 
Thus, /ii > or /i2 > can be assumed. Therefore, case 
III can be ignored without loss of optimality. 

. Case IV: f x = and f 2 = 0, then Ai > and 
A2 > 0. In this case, 9 belongs to the solution set 
for the equations given in ( |20] > and ( |2T| , which are 
obtained by satisfying the constraints in ( fTo} and ( fTTj ) 
with equality. ( |20| ) and ( |2Tj i correspond to |ATi| = C£ 
and \X 2 \ = C* 2 , respectively, which can be assumed 
without loss of generality due to Proposition [T] 
This completes the proof of Theorem [3] 

Appendix B 
Solving 9^ and b in Theorem[3] 

In this appendix, we show the steps to solve 6^. 0g can be 
obtained similarly. The unknown variables in ( [20] > are 8 and 
t. After some manipulations, ( f2"0"] > can be written as 

tcos-q + dit 2 + d 2 =0 (25) 
t cos(?7 + ui) + d 3 t 2 + d 4 = (26) 

where 

U = 2(022 + 011 - 012 " 02l), V = + 2(012 ~ 01l) (27) 
, A «l!^12| 4 -l ^ Qi|/lii| 4 C^+6i 

1 2a 1 |/i 11 | 2 IM 2 Q 1 ' 2 2a 1 |/ ll i| 2 |/n 2 | 2 Q 1 ' 

, A l^ 2 | 2 = ( Q 2|^2l| 4 -l)Q 2 +b2 

3 21^1^' 4 2a 2 |/ i2 i| 2 |/ l 22| 2 Ci 1 ' 
From |26|, we have 

t sin 77 sin w = t cos 77 cos w + d^t 2 + di =>■ 

t 2 (l — cos 2 77) sin 2 a; = (icos7icosu; + d 3 t 2 + d^) 2 (28) 

Solving cos rj from ( p5| >, we have 

cos 77 = — (dit 2 + d 2 )/t. (29) 



Substituting ( |29] i into ( |28| i gives the following fourth order 
polynomial equation with respect to t: 

[t 2 — {d\t 2 + d 2 ) 2 ] sin 2 to — [(d 3 — d\ cosuj)t 2 + d 4 — d 2 cosw] 2 

Since the above equation only has t 2 terms, it can be trans- 
formed to the following quadratic equation by setting z = t 2 , 

eiz 2 + e 2 z + e 3 = 0, (30) 
where ei =d\ + d\ — 2d\d 3 cosw, 

e 2 =2(«fiG?2 + d 3 di) — 2(dic?4 + ^2^3) cos a; — sin 2 ui, 
e 3 =d\ + d 2 — 2d 2 d^ cos ui. 

Then z can be easily solved. Since z — t 2 and t is the 
magnitude of C X2 , only the solutions of z that are real 
and satisfy < z < C* 2 need to be kept, whereby the 
values for t are obtained. For those values of t satisfying 
|(g?i£ 2 + d 2 )/t\ < 1, we can get the value for 77 based on ( f2"9] >, 
i.e., 77 = arccos[— (dit 2 +d 2 )/t] or 77 = 27r — arccos[— (dit 2 + 
d 2 )/t). Then 9 can be obtained from |27|. If no such solutions 
exist, then 9^ is set to empty. 
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